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Two kinds of maps that describe evolution of states of a subsystem coming from 
dynamics described by a unitary operator for a larger system, maps dehned for fixed 
mean values and maps defined for fixed correlations, are found to be quite different 
for the same unitary dynamics in the same situation in the larger system. An affine 
form is used for both kinds of maps to find necessary and sufficient conditions for 
inverse maps. All the different maps with the same homogeneous part in their 
affine forms have inverses if and only if the homogeneous part does. Some of these 
maps are completely positive; others are not, but the homogeneous part is always 
completely positive. The conditions for an inverse are the same for maps that are not 
completely positive as for maps that are. For maps defined for fixed mean values, 
the homogeneous part depends only on the unitary operator for the dynamics of 
the larger system, not on any state or mean values or correlations. Necessary and 
sufficient conditions for an inverse are stated several different ways: in terms of the 
maps of matrices, basis matrices, density matrices, or mean values. The inverse 
maps are generally not tied to the dynamics the way the maps forward are. A trace- 
preserving completely positive map that is unital can not have an inverse that is 
obtained from any dynamics described by any unitary operator for any states of a 
larger system. 

PACS numbers: 03.65.-w, 03.65.Yz, 03.67.-a, 03.67.Pp 

Keywords: open quantum dynamics, quantum error correction, reversible operation, inverse 
map, not completely positive map, affine map 



2 



I. INTRODUCTION 

How extensively can open quantum dynamics be described by changes of density matrices 
made by linear maps? How much do properties of the maps reflect physical properties of the 
dynamics? Does the dynamics determine the form of the map? Does the dynamics determine 
whether the map has an inverse? Does an inverse map correspond to the reversed dynamics? 
Does it correspond to any dynamics at all? When we go from a map to its inverse are we 
still in the physical realm of operations that can be done to the system through dynamical 
interactions with other systems? Or are we only in the mathematical realm of maps? Can 
these questions be handled by using only completely positive maps? 

The most basic questions are about the connection between maps and dynamics. They 
are brought out here by the fact that the same open dynamics can be described by two 
different kinds of maps. These questions are put in sharper focus when inverse maps are 
considered. 

Quantum information theory has given emphasis to completely positive maps, and quan- 
tum error correction works with partial inverses that also are completely positive maps 



ISj. Completely positive maps describe evolution of states of a system that in- 
teracts with another system but is initially not correlated with the other system. These 
maps can be engineered simply by bringing the uncorrelated systems together and letting 
them interact. More generally, evolution of states of a system in open quantum dynamics, 
a system that interacts with another system and initially may be correlated with the other 
system, can be described by maps that are not completely positive 

particular, in quantum information processing, unwanted correlations will bring in maps 
that are not completely positive. Moreover, the inverse of a map, indeed the inverse of a 
completely positive map, generally is not completely positive. 

Maps that are not completely positive can be used to learn things whose truth does not 
depend on the maps. For example, they have been used to describe how Lorentz transfor- 
mations of spin depend on momentum for a particle with spin 1/2 and positive mass, find 
that every Lorentz transformation completely removes the information from a number of 
spin density matrices 12|], find that finite Lorentz transformations can produce maximal 
entanglement of both the spins and the momenta from separable states of two particles with 
finite momenta [iJ], and find that the entanglement of two qubits can be increased from 



zero to maximal by an interaction on just one of them [l^. When the maps used as tools 
are set aside, the results about Lorentz transformations or open djTiamics remain. 

Things can happen in open quantum dynamics that are not described by completely 
positive maps. On the other hand, we will see here that working with maps that are not 
completely positive can bring in maps and properties of maps that are not connected with 
dynamics. They may not describe anything that can happen physically. No results will 
remain when the maps are set aside. 

Here we consider maps, completely positive or not, that describe evolution of states of 
a subsystem coming from dynamics described by a unitary operator for a larger system. 
Mathematically, this includes all trace-preserving completely positive maps; they all can 
be associated with unitary operators for larger systems ^. We consider two more general 
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and 



kinds of maps: those defined using only fixed mean values as map parameters 
those defined using correlations [9\. We carefully describe and compare them. Each gives a 
correct description of the change in time of every state in its domain. The domain can be 
different for the two kinds of maps. We will see from examples that the forms of the two 
kinds of maps can be quite different in the same situation for the same unitary dynamics in 
the larger system. The conditions for a map to have an inverse are very similar for the two 
kinds of maps. 

Maps that are not completely positive have limited domains. Does this make the condi- 
tions for an inverse more complicated for these maps than for completely positive maps? We 
look at each map in an affine form [3] where the homogeneous part is a completely positive 
linear map. All the different maps with the same homogeneous part have inverses if and only 
if the homogeneous part does. The conditions for an inverse are the same for maps that are 
not completely positive as for maps that are. For maps for fixed mean values, whether there 
is an inverse depends only on the unitary operator for the dynamics of the larger system; 
it does not depend on any state, mean values, correlations, or absence of correlations. For 
maps for fixed correlations, whether there is an inverse depends on the unitary operator for 
the dynamics of the larger system and on the state of the other part of the larger system; it 



does not depend on any correlations or absence of corre 
We will review, extend, and apply a framework 10 



ations between the two parts. 



ll| for describing open quantum 



dynamics as a map of density matrices that extends to a linear map of matrices, as a 
map of basis matrices, or as a map of mean values. Describing the map these different 



4 



ways gives alternative statements of necessary and sufficient conditions for an inverse: that 
the homogeneous part of the affine map does not map any nonzero matrix to zero; that 
hnearly independent basis matrices are mapped to hnearly independent matrices; or that 
the homogeneous part maps the hnear space of all matrices for the open system one to one 
onto itself, not into a smaller subspace. These conditions apply to both kinds of maps. For 
maps for fixed mean values, a necessary and sufficient condition for an inverse is that a 
linear map in a space of mean- value vectors that corresponds to the homogeneous part does 
not map any nonzero mean- value vector to zero, or that a multiple of the identity matrix 
is the only density matrix that the homogeneous part maps to a multiple of the identity 
matrix. For a qubit, this is the condition that the center is the only point in the Bloch 
sphere that is mapped to the center. For maps for fixed correlations, the condition is that 
the zero mean-value vector and another vector are not both mapped to the same vector, 
or that a multiple of the identity matrix and another density matrix are not both mapped 
to the same density matrix. For a qubit, this is the condition that the center and another 
point in the Bloch sphere are not both mapped to the same point. 

For maps for fixed mean values, a necessary and sufficient condition for an inverse can be 
stated simply in terms of the unitary operator for the dynamics in the larger system: that 
it does not map all the matrices for the subsystem into a subspace of the linear space of 
matrices for the larger system that is linearly independent of a subspace of matrices for the 
subsystem. For maps for fixed correlations, an example will show that a condition for an 
inverse can not generally be stated in terms of the dynamics alone. 

We consider an example of the open dynamics of one qubit that interacts with another. 
For both kinds of maps, we find that the homogeneous part is a completely positive map 
with an inverse that is not completely positive. The inverse map is not tied to the dynamics 
the way the map forward is. We look at the reason for this most closely for the map for 
fixed mean values. In that case, the evolution of states of the one qubit that comes from the 
dynamics of the two qubits going forward in time is described by a single map specified by 
a fixed set of mean values of quantities that involve the other qubit. The reversed dynamics 
of the two qubits gives many separate maps that take the states of the one qubit back in 
time. These maps generally are different for different states of the one qubit because they 
depend on mean values involving the other qubit that are changed by the dynamics so that 
at the later time they are different for different states of the one qubit. These separate maps 
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do not join to form a single inverse map. The inverse map is not obtained this way. It does 
not come from the reversed dynamics the way the forward map comes from the forward 
dynamics. There is a similar disconnect of the reversed dynamics and the inverse map for 
fixed correlations. 

These are examples of a more general result. We show that a trace-preserving completely 
positive map that is unital can not have an inverse that is obtained from any dynamics 
described by any unitary operator for any states of a larger system. 

Our results provide answers to the basic questions asked in the first paragraph. Open 
quantum dynamics generally can be described by changes of density matrices made by linear 
maps. We describe two different kinds of maps that can be used. Properties of the two kinds 
of maps can be quite different for the same dynamics. The form of the map as a whole is 
different. For one kind of map, the dynamics determines whether the map has an inverse; 
for the other kind it does not. An inverse map generally does not correspond to the reversed 
dynamics. It might not correspond to any dynamics at all. It does not, for example, for any 
trace-preserving completely positive map that is unital. Inverses generally exist only in the 
mathematical realm of maps; the connection to the physical realm is broken when inverses 
are taken. 
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II. FRAMEWORK 

We use a framework that has been nearly completely developed already 1^, 
Consider two interacting quantum systems S and R, both described by finite matrices: 
NxN matrices for S and MxM for R. We use the basis matrices F^q, Fq^, and F^^ described 
previously [ll|. The F^o for = 0, 1, . . . N'^ — 1 are A^^ Hermitian matrices for S such 
that Fqq is Is, the unit matrix for S, and 



Trs [F^oF, 



N6 



(2.1) 



This implies that the F^o are linearly independent, so every matrix for S* is a linear com- 
bination of the F^Q. For example, the F^q for ^ = 1,2,... A^^ — 1 could be obtained by 
normalizing standard generators j.l6|] of SU{N). The Fq,^ for z/ = 0, 1, . . . — 1 are 
Hermitian matrices for R such that Fqq is Ir, the unit matrix for R, and 



Tr^ [Fo^Fo,] = M6^ 



(2.2) 
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Every matrix for i? is a linear combination of the Fq^. We use notation that identifies F^q 
with F^o ® Ifl and Fq^ with I5 (8> Fq^, and let 

= F^o ® Fo,. (2.3) 

Every matrix for the system of 5* and R combined is a linear combination of the F^,^. 

We follow common physics practice and write a product of operators for separate systems, 
for example a product of Pauli matrices S and H for two different qubits, simply as SS, not 
S (g) S. Occasionally we insert a (8> for emphasis or clarity. 

The matrices F^o for positive /i and Fqi, for positive u are generalizations of Pauli matrices 
(and like the Pauli matrices they have zero trace). We use them to describe density matrices 
the way we use Pauli matrices to describe density matrices for qubits. If 11 is a density matrix 
for the system of S and R combined, then 

n = + E {Fc.o)F^o + E E {F,u)F,,] (2.4) 

^^^^^ 0=1 fj.=0 u=l 

and the density matrix p for 5* is 

p = Tr^n = -[Is + E (F,o)F,o] (2.5) 

a=l 



SO that 



and in particular 



= Tr [F,^n] (2.6) 



(F,o) = Tr5 [F,o Tr^n] = Tig [F^op] • (2.7) 
If f/ is a unitary matrix, then 



with the tap-^v elements of a real orthogonal matrix lOj, so that t^^.^^ is t^y^ap- Since 
;7U[/ and UlU^ are 1, 

^00; 0/3 = ^Oa 1^0/3 ; ; 00 = ^«0 '5/30- (2-9) 
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III. DYNAMICS 



We consider evolution described by a unitary matrix U for the system of 5* and R com- 
bined. In the Heisenberg picture, each matrix A that represents a physical quantity for S 
is changed to U'^AU. Its mean value is changed to (f/Mf/). These changes of mean values 
determine the change of the state of S. In particular, the state of S is described by the 
mean values (-Fao) for positive a, which determine the density matrix p in Eq. fl2.5p . Taking 
mean values in Eq. (12.81) gives 

) (3.1) 

11=1 11=0 u=l 

for a = 1, 2, . . . A^^ — 1. This is for one state of S that is part of a state of S and R combined. 



IV. MAPS FOR FIXED MEAN VALUES 



We would like to describe the open dynamics of S with a single map that applies to a set 
of different states of S. We will see two ways this can be done. In Section VII, we consider 
maps for fixed correlations [9)]. Here we describe maps for fixed mean values [lo|,llll- 

We define a map Vt of mean values (-Fqo) for positive a by letting 



fi((F,o)) = {U^F^oU) (4.1) 

with the {F^y) for positive z/ in Eq.f l3.1l) held fixed. We consider the {F^y) for positive v that 
participate in Eq. fl3.1l) to be map parameters. They describe the effect of the dynamics of 
the larger system of 5* and R combined that drives the evolution of S. There may be {F^y) 
for positive v that are not map parameters; they have no effect on the change of states of 
S because in Eq.f l3.ip they are with t^Q.^y that are zero. Which {F^^) are map parameters 
depends on U . 

The map VL applies to all the states of S described by mean values (-Fqo) that are com- 
patible with the fixed map parameters {F^y) in describing a possible initial state of S and R 



combined. This set of states of S is called the compatibility domain of the map [10|, lUl] . For 
each state of 5* in the compatibility domain, the state of S and R is described by a density 
matrix IIq given by Eq. fl2.4p for the (-F^o) that describe the state of S and the fixed map 
parameters (F^^). 
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The mean value (A) = TTs[Ap] for a matrix A for S is changed to 

(U^AU) = Tr[f/Mf/nn] = Tts[A TTR[UUnU% (4.2) 
so the Schrodinger picture for S is that the density matrix p for S is changed to^ 

n{p) = TiRpUnU^] = L{p) + K (4.3) 

where 

M 



L{Q) = T,n[UQ^U^] (4.4) 



-n 1 Af2-lAf2_l 

K = TmpiUn - p -|)f/t] = E E {F,.)Tr r[UF,M']. (4.5) 

^'^ ^^^'^ ^1=0 u=l 

The L part is a completely positive linear map that applies to any matrix Q for S, density 
matrix or not. It has the property that L{ts) is Is- The map L depends on U but does not 
depend on the state of R or on mean values of any quantities that involve R. The matrix 
K is the only part of Q{p) that depends on the state of R or on mean values that involve 
R. It does not depend on the state of S. It depends on U and on the map parameters {F^j^^) 
for positive z/, but not on the (Fao). With the same K, the Eg. (14. 3 p defines a map Q that 
applies to different density matrices p representing different states of S. 

Different map parameters (Ff^u) specify different maps. Each map Q applies to different 
states of 5* described by different (-Fqo)- For each map Q there is one NxN matrix K. 

The map Q of density matrices p, described by Eqs.f l4.3l) - (l4.5p . extends to a linear map 
of all matrices for S. The map of the basis matrices F^o to 

Q{ts) = ts + NK, Q{F^o) = L(F„o) (4.6) 

for positive a is equivalent to Eqs. fl4.3l) -f l4.5p for a set of p described by Eq. fl2.5p with variable 
(Fao). Since K is the same for all p, it cannot come from the terms with variable coefficients 
{Fao)', it can only be part of ^2(15). For each matrix Q for S, density matrix or not, the 
map is that 

Q{Q) = L{Q) + KTts{Q). (4.7) 

The change of states of S can be described equally well by the map Cl of the mean values 
(Fao) for positive a. For each state, these (Fao) are the components of a vector in a space 



^ The notation is that for each map X of matrices, the corresponding map of mean values is X. 
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of A^^ — 1 dimensions. Consider the change of states of 5* described by L. Because -^(15) is 
I5, the map L for this change of states extends to a hnear map of the space of mean-value 
vectors. To see this, we work out the equations that are needed. 

If mean values (-Foo)i are for a density matrix pi and (-Fao)2 are for a density matrix p2 
in Eq.f l2.5p . then the mean values for the density matrix 

p = ppi + (l-p)p2, (4.8) 

for p between and 1 are 

(F„o) = p(F„o)i + (1 -p)(F„o)2 (4.9) 

and, because 

L(p)=pL(pi) + (l-p)L(p2), (4.10) 

the mean values for L[p) are^ 

L((F„o)) =pL((F,o)i) + (1 -p)L((F,o)2). (4.11) 

In particular, when {Fao)2 is zero for all nonzero a, so 

{F^o) = p{Fao)i: (4.12) 

then P2 is U^/iV and L{p2) is Is/A^, because L{ls) is Is, so I/((Fao)2) is zero, and 

L{p{Fao)i) = L{{Fao)) = pL{{Fao)i) (4.13) 

for all nonzero a. When {Fao)2 is — (-^ao)i and p is 1/2, then {Fao) is zero, so p is ts/N and 
L(p) is ts/N, because L{ls) is I5, so L{{Fao)) is zero and 

L{-{F^o)i) = L{{F^o)2) = -L{{F^o)i) (4.14) 

for all nonzero a. These equations are all we need to see that L extends to a linear map of 
the space of mean- value vectors. 

V. INVERSE MAPS FOR FIXED MEAN VALUES 

We consider the open quantum dynamics of S described by a map Q using Eqs. fl4.3l) - fl4.5p 
and (14.71) for fixed mean values. For each L, the different maps Q for different K all have 
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inverses if and only if L has an inverse. This is immediately evident from the observation 
that they all are one to one if and only if L is. The inverse Q"^ of Q is that for any matrix 
Q for S 

n-\n{Q)) = L-\n{Q) - kttsMQ)]) (5.1) 

because Tr5[f2((5)] = Tr5'[(5]. The inverse maps generally take density matrices to density 
matrices only in limited domains for limited sets of states of S. In this they are like the 
maps for the open quantum dynamics going forward 10|, lUj . They apply where they are 
meant to be used. 

These statements are for all the maps Q for fixed mean values. These maps are made to 
describe evolution of states of S coming from dynamics described by a unitary operator U 
for a larger system of 5* and R combined. Some of these maps Q are completely positive; 
others are not, but L always is completely positive. The conditions for an inverse are the 
same for maps that are not completely positive as for maps that are. Whether there is an 
inverse depends only on L, which depends only on U. It does not depend on the state of R 
or on any mean values that involve R. 



A necessary and sufficient condition for an inverse 17| is that L does not map any nonzero 



matrix to zero. We can see from Eq.f l4.4p that L{Q^) is L{Qy for any matrix Q for S. If 
L[Q) is zero, so is ^(Q"''). Hence the necessary and sufficient condition for an inverse is that 
L does not map any nonzero Hermitian matrix to zero. We can see from Eq. (l4.4p how this 
depends on U. 

I — I 

A necessary and sufficient condition for an inverse [17] is that L maps the basis matrices 
Fao for S to linearly independent matrices. In fact, the necessary and sufficient condition 
for an inverse is just that the matrices L{Fao) for positive a are linearly independent, be- 
cause L{ls) is linearly independent of them, because TTs[L{Fao)] is zero for positive a but 
Tr5[L(l5)] is not zero, because L does not change the trace of F^o, which we can see from 
Eq. (14.41) . For each L, the matrices L{Fao) for positive a are the matrices fl{Fao) of Eqs. (14.61) : 
they are the same for all the different maps Q for different K. 

Every matrix for S* is a linear combination of the Fao- A necessary and sufficient condition 
for an inverse is that L maps the linear space of all matrices for S one to one onto itself. The 
alternative, the necessary and sufficient condition for no inverse, is that L maps the Fao to 
linearly dependent matrices and maps the linear space of all matrices for S into a subspace 
of itself that has fewer dimensions. Then there is a nonzero subspace of the linear space of 
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all matrices for S that is linearly independent of the subspace of all L{Q) for matrices Q 
for S. From Eg. (14.41) . we can see that this gives a necessary and sufficient condition for an 
inverse in terms of U : that in the linear space of all matrices for 5* and R combined there 
is no nonzero subspace of matrices for S that is linearly independent of the subspace of all 
UQU^ for matrices Q for 5*. 

In terms of L, shown in Section IV to be a linear map of the space of mean- value vectors, 



a necessary and sufficient condition for an inverse [17|] is that L does not map any nonzero 
mean-value vector to zero. For a qubit, this is the condition that the center is the only 
point in the Bloch sphere that L maps to the center. Since the components of a mean- value 
vector are the mean values (-Fqo) for positive a, the necessary and sufficient condition for an 
inverse is that the {Fao) for positive a are not all zero after the dynamics if they were not 
all zero before. From Eg. (12.51) . we see that in terms of density matrices the necessary and 
sufficient condition for an inverse is that L{p) is Is/A^ only if p is l^/A^. 

VI. EXAMPLES FOR FIXED MEAN VALUES 

We consider two qubits, described by Pauli matrices Si, S2, S3 for S and Si, S2, H3 
for R, so Fjo is Sj and F^k is S^, which implies Fjk is SjSfc, for j, k = 1,2, 3. The density 
matrix for the two qubits is 

\ J=l k=l j,k=l J 

and the density matrix for 5* is 

p = Tr«n = i + I . (6.2) 



Let 



We can easily compute 



t/^Sif/ = Sie^*^^^'^^ 

= Si C0S7 — S2S3 sin7 (6.4) 
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using the algebra of Pauli matrices, and similarly 

f/^S2f/ = S2 COS7 + S1H3 sin7 

f/l'Sgt/ = S3. (6.5) 

Interchanging U and W has the same effect as changing the sign of 7. Thus we see that 

L{ls) = Is 
L(Si) = S1COS7 
L(S2) = S2COS7 

HJ^s) = S3. (6.6) 
The map L is specified by these Eqs. fl6.6p . or by the corresponding map of mean values 

L((Si)) = (Si)cos7 
L((S2)) = (S2)cos7 

L((E3)) = (S3), (6.7) 

or by 

L{Q) = cos(7/2) Q cos(7/2) + S3 sin(7/2) Q S3 sin(7/2) (6.8) 

for any matrix Q for S, which gives Eqs. fl6.6p . 

The necessary and sufficient condition for L to have an inverse is that cos 7 is not zero. 
Then the inverse map L^^ is specified by 

L'\ls) = Is 

L^"^(Si) = Si/ C0S7 

L"^(S2) = S2/ C0S7 

L-i(S3) = S3, (6.9) 

L-i((Si)) = (Si)/cos7 
L-i((S2)) = (S2)/cos7 

L-i((S3)) = (S3), (6.10) 



or 



T-i(n\ ^ cos(7/2) cos(7/2) sin(7/2) sin(7/2) 

C0S7I w|cos7| w|cos7| -i/|cos7| 
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with the upper signs for positive cos 7 and the lower signs for negative cos 7. Clearly is 
not completely positive. It does not map all positive matrices to positive matrices. 

Here we considered the map fl that is just L, with no K. From Eqs. fl6.4p and (16.51) . we 
can see that the map Cl of mean values generally is 

r2((Si)) = (Si) C0S7 - (S2S3) sin7 
Cl{{T.2)) = (S2) C0S7 + (S1H3) sin7 

fi((S3)) = (S3). (6.12) 

It depends on (E2S3) and (S1H3) as well as on U. The map Q is just L, with no K, when 
(E2S3) and (S1S3) are zero. In general, we need 

2K = [-(E2S3)Ei + (SiS3)S2] sin7 (6.13) 
to get the mean values in Eqs. (l6.12l) with 

L{p) + K = Q{p) = l[n{is) + j2{j:,)L{j:,)]. (6.i4) 

VII. MAPS FOR FIXED CORRELATIONS 

Now we describe maps for fixed correlations j^. We return to Eq. fl3.ip . which gives the 
changes of the mean values (-Fqo) for one state of S. From there we take a different route to 
describe the open quantum dynamics of S with a single map that applies to a set of states of 
S. The result is a different map, sometimes an infinite number of different maps, replacing 
the map Q that we obtained for fixed mean values. 

We use the correlations 

r,u = {F^u)-{F^o){Fou) (7.1) 

for positive /i, u. We define a map $ of the mean values {Fao) for positive a by using 
Eq.f l7.1l) to substitute for {Ffj^i,) in Eq. fl3.ip and letting 

<l((F„o)) = {U^F^oU) 

= X! taO;^io{F^o) + ^ taO;Ou{Fou) 
fj,=l v=l 

N'^-lhP-l Af2-1A/2_1 

(7.2) 

11=1 v=l ^=1 v=l 
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for fixed mean values (Fq^) or density matrix 

PR = j^[tR+ E {Fou)Fo.] (7.3) 

for the state of R, and fixed correlations T^i,. We consider the (Fq^) and T^u that participate 
in Eq. fl7.2p to be map parameters. They describe the effect of the dynamics of the larger 
system of S and R combined that drives the evolution of S. There may be (Fq^) and F^j/ 
that are not map parameters; they have no effect on the change of states of S because 
in Eg. (17.21) they are with tao-ou or tao-fiu that are zero. Which {Fq^) and V^^ are map 
parameters depends on U . 

This map $ applies to all the states of S described by mean values (-Fqo) that are com- 
patible with the fixed map parameters {Fq^) and F^,^ in describing a possible initial state 
of S and R combined. We call this set of states of 5* the compatibility domain of the map 
$. It may be smaller than the compatibility domain of a corresponding map for fixed 
mean values, because for $ the {Fao) have to be compatible with a larger number of fixed 
parameters. Generally there is an added parameter for $ when there is a /i and a v for 
which tao-ou is zero for all a but tao-^u is not. Then {Ff^u) is a parameter for Q but {Fq^,) 
is not, and both F^j^ and (Fq,,) are parameters for $ in all but the exceptional cases where 
there is a cancellation of tao -^u{Fou) with tao-^o- In Sections VI and IX we will see examples 
where Q depends on two parameters and $ depends on three. 

For each state of S in the compatibility domain, the state of S and R is described by a 
density matrix 

n$ = T77j[l+ E {Fao)F.o+ E iPou)Fo. 

Af2-1 A/2-1 Ar2-l A/2-1 

+ E E {F,,){F,,)F,,+ E E (7.4) 

obtained by using Eq. (l7.ip to substitute for {F^y) in Eg. (12.40 . 

The mean value {A) = TTs[Ap] for a matrix A for 5* is changed to 

{U^AU) = Tip^AUU^] = Tts[A Ttr[UU^U^] (7.5) 

so now the Schrodinger picture for S is that the density matrix p for S is changed to^ 



$(p) = TvrIUU^W] = D{p) + C 



(7.6) 
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where 

D{Q) = 1:ir[UQprU^] (7.7) 
C = TinMn^ - ppn)U^] = -— J2 r,,TTn[UF,M^]. (7.8) 

^^-''^ fj,= l u = l 

The D part is a completely positive linear map that applies to any matrix Q for density 
matrix or not. The map D depends on U and on the state of R represented by the density 
matrix p^. It does not depend on the correlations T^i, between S and R. The matrix C is 
the only part of $(p) that depends on the correlations T ^y. It also depends on U , but it 
does not depend on the state of S described by the mean values (F^o) or density matrix p 
or on the state of R described by the mean values {Fq^) or density matrix pR. 

Now the map parameters are the {Fq^) and T ^y. Different {FqJ) and r^^, specify different 
maps. Each map $ applies to different states of S described by different (-Fqo)- For each 
map $ there is one density matrix p^ and one NxN matrix C. 

The map $ of density matrices p described by Eqs.( l7.6l) -( l7.8p extends to a linear map of 
all matrices for S. The map of the basis matrices F^o to 

^Is) = ls + NC, ^F^o) = D{F^o) (7.9) 

for positive a is equivalent to Eqs.f l7.6l) -f l7.8p for a set of p described by Eq.f l2.5p with variable 
(Fao)- Since C is the same for all p, it cannot come from the terms with variable coefficients 
(Fao); it can only be part of $(1^). For each matrix Q for S, density matrix or not, the 
map is that 

^Q) = DiQ) + CTTsiQ). (7.10) 

Let E be the map defined by 

E{Q) = D{Q) - D{ls)'lTs{Q)lN + ts'lTs{Q)/N (7.11) 

for all matrices Q for S. Then E is & linear map, and E{'ts) is l^, so the corresponding 
map E of mean values is a linear map of the space of mean-value vectors, as was shown for 
L in Section IV. 

The maps $ include the familiar completely positive maps that have long been used 
for initial states with no correlations. When the correlations T are zero, C is zero and 
Eqs. ( 17.6(1 and (17.71) are the familiar equations that define the familiar completely positive 
maps (5|. When maps defined with fixed mean values are used, the familiar completely 
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positive maps are separate different maps, alternative maps that are options when there are 
no correlations in the map parameters |10| . 



VIII. INVERSE MAPS FOR FIXED CORRELATIONS 

Inverses of maps for fixed correlations parallel those for fixed mean values. For each D, 
the different maps $ for different C all have inverses if and only if D has an inverse. The 
inverses are described by Eq. fl5.ip with fi, L and K replaced by $, D and C . Again, the 
inverse maps generally take density matrices to density matrices only in limited domains for 
limited sets of states of S. In this they are like the maps for the open quantum dynamics 
going forward. They apply where they are meant to be used. 

These statements are for all the maps $ for fixed correlations. These maps are made to 
describe evolution of states of 5* coming from dynamics described by a unitary operator U 
for a larger system of S and R combined. Some of these maps $ are completely positive; 
others are not, but D always is completely positive. The conditions for an inverse are the 
same for maps that are not completely positive as for maps that are. Whether there is an 
inverse depends only on D, which depends only on U and on the state of R described by 
the density matrix p^. It does not depend on the correlations, or absence of correlations, 
between S and R. 

A necessary and sufficient condition for an inverse is that D does not map any nonzero 
Hermitian matrix to zero, or that the matrices D(Fq,o) for positive a are linearly independent. 
For each D, the D{Fao) for positive a are the $(Fao) of Eqs. (17.91) : they are the same for 
all the different maps $ for different C . The discussion for fixed mean values in Section V 
applies almost unchanged with L replaced by D. One change is that -D(ls) is generally not 
I5. In terms of the map of mean-value vectors, a necessary and sufficient condition for an 
inverse is that E does not map any nonzero mean-value vector to zero, which means that D 
does not map the zero mean-value vector and another mean-value vector both to the same 
vector, which means that D does not map t-s/N and another density matrix both to the 
same density matrix.^ For a qubit, this is the condition that D does not map the center and 
another point in the Bloch sphere both to the same point. For maps for fixed correlations, 
there is no necessary and sufficient condition for an inverse in terms of U alone, as there is 
for maps with fixed mean values; we will see an example where the condition for an inverse 
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involves both U and the state of R. 



IX. EXAMPLES FOR FIXED CORRELATIONS 

Again, we consider two qubits, described by Pauli matrices Si, S2, S3 and Hi, S2, H3, 
and an open quantum dynamics of the S qubit that comes from the dynamics described by 
the unitary operator U of Eq. fl6.3p for the two qubits. Now we let each map be specified by 
fixed mean values (S^) and fixed correlations 

r,fc = (s,s,)-(s,)(s,) (9.1) 

for j, k = 1, 2, 3. Taking mean values in Eqs. (16.41) and (16. 5p and using Eqs. (l9.ip to substitute 
for (S2S3) and (S1S3), gives the map of mean values 

<l((Si)) = (Si) C0S7 - (S2)(S3) sin7 - r23sin7 
<I>((S2)) = (S2) C0S7 + (Si)(H3) sin7 + ri3sin7 
<I>((S3)) = (S3). (9.2) 

From calculations like those in Eqs. (l6.4p and (16. Sp . with the sign of 7 changed to interchange 
U and W , we find that 

D{ls) = Is 

D(Si) = Si cos 7 + S2 (S3) sin 7 

D{T.2) = S2 C0S7 - Si(E;3) sin7 

Z}(S3) = S3 (9.3) 

C = ^(ri3S2-r23Si)sin7. (9.4) 

We can see that these give a map $ that gives the same changes of density matrices p as 
the map $ of mean values described by Eqs. (l9.2p . The map $ is completely positive only if 
Fis and r23 are both zero; only then is the matrix 

$(i[]l5 + S3]) = i[l5 + S3 + ri3S2sin7-r23Sisin7] (9.5) 

a positive matrix. 

The condition for an inverse is that the determinant 

det = (cos7)2 + (S3)2(sin7)2 (9.6) 
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is not zero, which means that cos 7 and (S3) are not both zero. If cos 7 is zero, the condition 
for an inverse is that (H3) is not zero. Whether there is an inverse depends on both U and 
the state of R. The inverse map of basis matrices is 

B-'(ls) = Is 

det det 

det det 
d-\j:s) = S3. (9.7) 

The corresponding map of mean values is 

D-Hi^s)) = (S3). (9.8) 



This gives 



[D-i((Si))]^ + [D-^m))r = ^^'^'^J^"^ (9-9) 



2 



which shows that is not completely positive unless (S3) is 1 or sin 7 is zero. From 
Eqs.f l9.2l) we see that if sin 7 is zero, then $ can either do nothing or just change the signs 
of (Si) and (S2). If (H3) is 1, then the state of i? is a pure state, so 11$ is a product pp^ 
for any state of S in the compatibility domain, which means there are no correlations, and 
the map described by $ or $ is just rotation by 7 around the z axis. 

X. DYNAMICS DISCONNECTION 

For the larger system of S and R combined, we can consider both the dynamics described 
by U that takes states forward for an interval of time, and the reverse dynamics described by 
W that takes states backwards through the same interval. We consider maps for constant 
mean values first. In our example for constant mean values, we assume now that (S2H3) 
and (S1S3) are zero at the start of the interval, so the map Q that describes the evolution of 
states of S going forward through the interval is just L, with no K. Since the change from U 
to W is just the change of sign of 7, the maps that describe the changes of states of S going 
backwards through the interval have the same L, described by Eqs. (l6.6p or (16.81) . and have K 
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matrices described by Eqs. fl6.13p . which give maps of mean values described by Eqs.f l6.12!) . 
with sin 7 changed to — sin 7 and (S2S3) and (S1S3) changed to the values (S2S3)' and 
(E1S3)' at the end of the interval, which, from calculations like those of Eqs. fl6.4p and (16.51) . 
can be seen to be 

(S2S3)' = (Si)sin7 

(S1S3)' = -(S2)sin7 (10.1) 

with (Si) and (S2) mean values at the start of the interval. After the changes forward and 
backward, the mean values (Sj) are returned, as we know they must be, to 

(51) " = L((Si)) C0S7 + (S2S3)'sin7 

= (Si)(cos7)' + (Si)(sin7)' = (Si) 

(52) " = L((S2)) C0S7 - (SiH3)'sin7 

= (S2)(cos7)' + (S2)(sin7)2 = (S2) 

(53) " = L((S3)) = (S3). (10.2) 

The maps that describe the changes of states of 5* going backwards in time are different for 
different values of (Si) and (S2); their K matrices are different because the (S2H3)' and 
(S1S3)' given by Eqs. fllO.ip are different. These maps generally are different for states of S 
that are different at the start of the interval. They do not comprise a single inverse map for 
all the states of S that are mapped forward. 

The inverse map in our example is not obtained from U'^ and equations like Eqs. (14.31) - 
(14. 5p . It is not tied to the dynamics the way L is. In fact, the map L^^ can not be obtained 
the way L is, as a map for fixed mean values, from any dynamics described by any unitary 
operator for any states of a larger system of S and R combined. There are two ways to see 
this. One is that if L^^ were obtained that way, it would have to be completely positive, 
because its K would be zero, because L^^(ls) is Is; but L^^ is not completely positive. 
Another way to see that L^^ is not obtained from dynamics, described next, will apply to 
our example for fixed correlations as well. 
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XI. NO DYNAMICS POSSIBLE 



There is no way that an inverse can be obtained from dynamics for any unital trace- 
preserving completely-positive map that describes truly open dynamics. A unital map is a 
map that takes the unit matrix to the unit matrix, so it maps the density matrix for the 
completely mixed state to itself. 

Theorem. If a completely positive map preserves the trace, is unital, and is not simply a 
unitary transformation for the subsystem 5* alone, then it does not have an inverse that is 
obtained from any dynamics described by any unitary operator for any larger system of S 
combined with another system R. 



Proof. A trace-preservmg completely pos.tive map ca. be obtained from dynamics as a map 

for fixed correlations that are zero [5^ . It is a map D described by Eq. (17. 7p for some R and 
U. If the map is unital, then D{ls) is ts, and the map E defined by Eq.f l7.11i) is D, so D 
is a linear map of the space of mean- value vectors.^ If the map has an inverse D~^, then 
D^^ is the inverse of D and is a linear map of the space of mean- value vectors. The density 
matrix p for a state of S, described by Eq. (12.51) . is changed by D to 



1 



a=l 



;ii.i) 



with D{{Fao)) = {WFaoU) calculated from Eq.i^. Since D is a linear map of the space 
of mean- value vectors, each D{{Fao)) given by Eq. (l3.1l) is a linear combination of the (F^q) 
alone. This implies that 



Tr 



iD{p)y 



< Tr 



P 



To see this, let 

= (Fao) for « = 1, 2, . . . - 1 
and Xa/3 = for other a, /3, and let 



(3 = 



A/2-1 

E E 

a=0 /3=0 



'11.2) 



11.3) 



;ii.4) 



Then D{{F^q)) is D{x^o) for positive /i, because t^o;00 is zero, and 



Af2-1 Af2-l A/2-1 

E (^((^.o))) < E E [Dix^u] 

fj.=l /i=0 i/=0 



Ar2-1 A/2-1 

E E iXap) 
a=0 p=0 



Af2-1 

E (^-o)^ 



;ii.5) 
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because the t^iy ■ ap are the elements of a real orthogonal matrix. This proves the inequality 
flll.2p . If D^^ comes from dynamics (with both R and U generally different than for D), 
we can see similarly that 



Tr 



P 



< Tr 



(Dip))' 



11.6) 



Pure states are mapped to pure states. For each vector \tp) of length one, there is a vector 



l^jj ) of length one such that 



Absolute values of inner products of vectors are not changed: if 

P=^IV^)(V^I + ^I0)(0| 



for two vectors and |0) of length one, then 



^ + ^l(^l0)r = Tr 



P 



Tr 



iDip)y 



1 1 



:ii.7) 



;ii.8) 



;ii.9) 



2 2' 

The trace of a product of Hermitian matrices is not changed: for one-dimensional projection 
operators \ip){ip\ and \4>){(p\ for vectors {ip) and |0) of length one. 



Tr 



1(^10)1^ = = Tr [D{\^){^\)D{\<i>)m 



11.10) 



and for any Hermitian operators F and G there are spectral decompositions 

^ = E/.#.)(^il' G = T.9km{<i^k\ (nil) 

j k 

with real numbers fj and and vectors \ifjj) and \(j)k) of length one, which give 

Tr[FG] = J2f,gkTTmm<l>k){A\] 
jk 

= Y.f,9kTrm^,){^,\)D{\<p,){<p,\)] = TT[D{F)DiG)]. (11.12) 

jk 

Equation (12.11) holds for the i5(F^o) as well as for the F^q. This implies that the A^^ matrices 
D(F^o) are linearly independent, so every matrix for S* is a linear combination of the D(-F^o)- 
The linear space of all matrices for S is mapped one-to-one onto itself. The set of all pure 
states is mapped one-to-one onto itself. Then Wigner's theorem says that the map is made 
by either a linear unitary operator or an antilinear antiunitary operator on the space of 
state vectors jisl. An antilinear antiunitary operator is not possible j2o|. The map is a 
unitary transformation for S alone. This completes the proof of the theorem. 
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In particular, in our examples, neither nor generally can be obtained from any 
dynamics described by any unitary operator for any states of a larger system of S combined 
with another system R. The only exceptions are when L is the identity map and when D is 
rotation of the S qubit around the z axis. 

XII. DISCUSSION 

Either way the map is defined, with only fixed mean values as map parameters or with 
fixed correlations, the map gives a correct description of the change in time of every state 
of S in its compatibility domain; the map describes the change that comes from the unitary 
dynamics in the situation described by the map parameters in the larger system of S and 
R combined. Our examples show that as a whole the map defined one way can be quite 
different from the map defined the other way. The map D described by Eqs.f l9.3l ) contains a 
rotation that is not in the map L described by Eqs.f l6.6( ) or (16. 7p . The compatibility domains 
of the two kinds of maps are different, and for states that are in both of their compatibility 
domains, maps Q and $ generally are different. In our examples, maps Q and $ give the 
same result only for the one set of values of (Si) and (S2) that make and f2((S2)) 

described by Eqs.f l6.12l) the same as $((^1)) and <I>((S2)) described by Eqs.f l9.2p : these are 
the (Si) and (S2) that are solutions of Eqs.( l9TT]) for the fixed parameters (S2H3) and (S1S3) 
of Q and the fixed parameters (S3), r23 and ri3 of $. 

Which kind of map is the better tool to use may depend on what the job to be done is. 
In particular, one compatibility domain may extend better than the other to include states 
of interest. The differences between the two kinds of maps do not reflect physical properties 
of the dynamics. Far from it; they are differences between two ways of describing the effect 
of the same dynamics of 5* and i? on a set of states of S, either a set of states compatible 
with fixed mean values or a set of states compatible with fixed correlations. 

Whether a map has an inverse is a property of the map that may be closely tied to the 
dynamics. The conditions for an inverse are very similar for the two kinds of maps. For 
maps defined with only fixed mean values as map parameters, we saw in Section V that a 
necessary and sufficient condition for an inverse can be stated just in terms of the unitary 
operator that describes the dynamics of S and R. For maps defined with fixed correlations, 
our examples in Section IX show that the dynamics alone does not always determine whether 
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a map has an inverse; it can depend on the state of R as welL The conditions for an inverse 
are the same for maps that are not completely positive as for maps that are. 

The dynamics is generally less closely tied to what the inverse is. We saw in Section X 
that when a map comes from dynamics, the inverse map generally does not come from the 
reversed dynamics. It might not come from any dynamics at all. We showed in Section XI 
that if a trace-preserving completely positive map is unital, it can not have an inverse that 
is obtained from any dynamics described by any unitary operator for any states of a larger 
system. When a map is an operation that can be done to the system, its inverse might not 
be. 

Many maps that are not completely positive describe operations that can be done to 
a system through unitary dynamical interactions with other systems. Many maps do not. 
Many inverses of maps that do, do not. All trace-preserving completely positive maps do 

We can appreciate how quantum error correction works with completely positive maps. 
Using only completely positive maps avoids maps that can not be realized physically. The 
recovery map can be completely positive because it is not a full inverse. It is an inverse only 
on a subspace of the code space [5] . If a trace-preserving completely positive map has a full 
inverse that is completely positive, the map comes from a unitary operator on the system 



alone, without any interaction with another system j4|. Completely positive maps do not 
describe everything that can happen in open quantum systems, but they are the maps that 
are most simply, completely, and certainly tied to dynamics. 
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